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Formulas and tables are given for the mean and standard deviation of ft = 20 logio r 
where the random variable r has the Nakagami-Rice distribution. This distribution is of 
interest in connection with the short-term fading characteristics of sonic received radio fields. 
A particularly simple formula for the mean of R is obtained in terms of the well-known 
exponential integral function — Ei(—x). Additional information concerning the median and 
interdecile range of R is also given. 

1. Introduction 

In certain radio propagation problems the field strength at a receiver may be approximated 
by (he vector sum of a constant vector and a Rayleigh-distributed vector [.Norton, Vogler, 
Mansfield, and Short, 1955). [f it is assumed that the Rayleigh-distributed vector has an rms 
amplitude k and that the constant vector lias an rms amplitude of unity, then the probability 
distribution of the resultant r of their sum as given by Nakagami [1940] and Rice [1944, 1945] 
may be put in the form 

P(r>r ) = 1 1"° r exp [ - (1 +/' 2 )/FJ/ (1 (2r/F>/r (1.1) 

where Iq(x) is the modified Bessel function of first kind and order zero. If the R ay leigh distrib- 
uted vector has an rms amplitude k h the constant vector has an rms amplitude k 2 , and the 
resultant of their sum is r h then (1.1) gives the distribution of r rJk- 2 as a function of k = k 1 /k 2 . 
The primary purpose of this paper is the determination of the mean R and standard devia- 
tion a R of R=20 log 10 v. Additional information concerning the median //((). 5) and interdecile 
range B(0.1)-i?(0.9) of R is taken from the paper by Norton et al. [1955]. The phase </> of the 
vector sum of a Rayleigh-distributed vector and a constant vector is discussed by Norton, 
Shultz, and Yarbrough [1952]. 

2. Calculation of R and <t r 

We will determine a R from the relation 

where R 2 is the mean of R 2 . Making use of (1.1), we find that 

!=p J" (20 log 10 r)r exp [-(l+r 2 )/k 2 }I Q (2r/k 2 )dr (2.1) 

7^=!J^ (20 log 10 rfr exp [-(l+r 2 )/k 2 }I (2r/k 2 )dr. (2.2) 



and 



In appendix 1 of this paper, these two integrals are evaluated and the three equations which 
follow are derived. 

#=(10 log 10 e)[-Ei(-l/k 2 )]. (2.3) 
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For k> 1, use 
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For &<1, use 

*%= (10 log w ^) 2 l-2^(-l/F)£*(l/F)+4j i 
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Id these equations, 
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re 



where 



-Ei(-x)= f 

J -a; e-^0 V J -a: J e / 



dt, x>0; 
<ft, x>0, 



with e>0; and Ci = 1.0990 19, a constant. The functions — Ei(—x) and E*(x) are well-known 
functions for which tables are available. See, for example, the book " Tables of Sine, Cosine, 
and Exponential Integrals/' U.S. National Bureau of Standards [1940]. The integrals 



f 



-Ei(-t) 
t 



dt 



and 



f 



j[-Ei(-t)]dt 



are evaluated in appendix 2 to six significant digits by numerical integration. 

Letting K=20 log 10 k, we can get the following asymptotic formulas from (2.3), (2.4), 

and (2.5). 

ForiT>20, i?^#+ 4iP-2.5068 (2.6) 
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-40 
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0.000 


-0.000 


0.061 


-35 


.001 


.000 


-.001 


.109 


-30 


.002 


.000 


-.002 


.194 


-25 


.007 


.000 


-.007 


.346 


-20 


.022 


.000 


-.022 


.616 


-18 


.034 


.000 


-.034 


.776 


-16 


.054 


.000 


-.054 


.980 


-14 


.086 


.000 


-.086 


1.238 


-12 


.136 


.000 


-.136 


1.569 


-10 


.214 


.000 


-.214 


1.999 


-8 


.335 


.001 


-.334 


2.565 


-6 


.524 


.017 


-.507 


3.279 


-4 


.813 


.107 


-.706 


4.036 


-2 


1.249 


.383 


-.866 


4.667 





1.894 


.953 


-.941 


5.094 


2 


2.808 


1.855 


-.953 


5.340 
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4.006 


3.064 


-.942 


5.465 


6 


5.448 


4.519 


-.929 


5.525 


8 


7.077 


6.155 


-.922 


5.551 


10 


8.835 


7.917 


-.918 


5.562 


12 


10. 679 


9.763 


-.916 


5.567 


14 


12. 580 


11.664 


-.916 


5.569 


16 


14.517 


13. 602 


-.915 


5.570 


18 


16. 477 


15. 562 


-.915 


5.570 


20 


18. 452 


17. 537 


-.915 


5.570 
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00 
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5.570 



430 



20.0 

10.0 
7.0 
5.0 







R(O.I 


-R(0.9) 
































































f °n y 




















































R(0.5) 


/ R 
























R(0.5) 


-R 








If J 












y / 












/ / 












/ I 












/ I 


















7 




/ 









0.2 7 — -<- 



-30 -20 -10 

Figure 1. K in decibels. 



and 



(7^4.3429 V1.6449-l/(2& 4 ). 

Fori£<-20, /7«4.;M29F6- 1/a2 

^«(hl418A:Vl+^72. 



(2.7) 
(2.8) 
(2.9) 



The magnitude of the error in these four expressions is less than 5(10) -4 . 

Values of the median, the mean, the difference between the mean and the median, and the 
standard deviation of R are given in table 1. These four quantities together with the inter- 
decile range of R are shown as functions of K in figure 1. The median and interdecile range are 
taken from the paper by Norton et al. [1955]. 

3. Appendix 1. Derivation of the Formulas for R and <r R 

If we change to natural logarithms and make the substitution r=ky in (2.1) and (2.2), 
we get 

i=2(20 log 10 e)e- 1/k2 [(\nk)A(k)+B(k)], 



and 
Here 



R 2 =2(20 log 10 e) 2 e- 1/ " 1 [(lnk) 2 A(k)+2(lnk)B(k) + C(k)}. 
A(k)= f ye-" 2 Io(2y/k)dy=[ P if a -'e^ 2 I (2y/k)dy'] , 

Jo LJo Ja=l 

*(*)=// (In y)ye-»%(2y/k)dy=^± f° tf°- 1 e->> 2 I (2y/k)dy'^ 



(3.1) 
(3.2) 
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and 

Using the Maclaurin series expansion for I (2y/k) and integrating term by term, we find that 

P e-»Y-%VylWy=^ iFx(a; 1; 1/F), 

JO 

where Y(x) is the gamma function and iF x (a; c; x) is a confluent hypergeometric function 
[Rainville, I960]. 

In order to complete the evaluation of A(k), B(k), and C(k), we make use of the known 
series expansion for —Ei(—x) and of the series transformation 

e~ x S a» T^S K ~P where b n =A n a (3.3) 

w =o ft! »=o ft! 

is the nth forward difference of a ; i.e., A°a =a , A 1 a = ai—a , etc. Omitting any further 
details, we find that R is given by (2.3), and that 

R 2 =(10 log 10 e) 2 /4(lni) 2 +4(lnJt)[ln (l/k 2 )-Ei(-l/k 2 )]+y 2 

+--2 T[7+ ln (l/F)-^(-l/F)]+2 g < ( „% .^ (gi)}- (3.4) 

Combining (2.3) and (3.4), we get (2.4) for a R . 

In order to prove (2.5), we derive a different expression for the series in (3.4). Observe 
that 



where 



(_l)«+i(a;)»+i / n n p. r . (_i)«+if» / . i\-i 
£i (n+l)!(n+l) V£ij7 Ji U=i (»+l)l Vi^i/J 

too / l\n+l r «+l / W 1\~| 
S ;,L , n ( 2 - ) =0.1827580. 
J=i (ft+l)!(ft+l) V^ j/J x =i 



Omitting any further details, except to note that use is made of the series expansions for 
—Ei{ — x) and E*(x), and that we again use the series transformation (3.3), we find that 

+ j; e JEi(-mt-2f i ^ d t+yln X+ ^f + 0. 
Using this expression together with (2.3) and (3.4), we arrive at (2.5). 

4. Appendix 2. Evaluation of the Integrals I ^~ — -dt and 

J x j[-Ei(-t)]dt 



The values of the integrals 



gd/k^j^^^dt and h(l/P)=f; /k ^l-Ei(- 



t)]dt, 



given in table 2 for various values of K, were obtained as follows. The ith entries of the 
table were computed in terms of the preceding entries by means of the formula 



j>><*-(X'+£; ■+...+£>)« 
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Table 2 



K 


Ilk - 


0(1/* 2 ) 


Ml//; 2) 


-0.96910013 


1. 2500000 


0. 0404747 


0. 123376 


-2 


1. 5848932 


. 0679414 


. 234919 


-4 


2. 5118864 


.0917103 


. 402303 


-(', 


3.9810716 


. 0971504 


.518481 


-6. 9897000 


5. 0000000 


. 0976706 


. 561387 


-8 


6. 3095733 


.0978116 


. 597259 


-9.0308999 


8. 0000000 


. 0978393 


. 626974 


-10 


10. 000000 


. 0978428 


. 649645 


-10.791812 


12. 000000 


. 0978432 


. 665007 


-12 


15. 848932 


. 0978432 


. 683943 


-14 


25.118864 


. 0978432 


. 706139 


-16 


39.810716 


. 0978432 


. 7203s:-; 


-18 


63. 095734 


. 0978432 


. 729470 


-20 


100. 00000 


. 0978432 


. 735245 


-25 


316. 22776 


. 0978432 


. 742039 


-30 


1000. 0000 


. 0978432 


.744190 


-35 


3162. 2776 


. 0978432 


. 744880 


-40 


10000. 000 


. 0978432 


. 745090 


— 00 


00 


. 0978432 


. 745190 



where each of the integrals on the right-hand side of the equation was evaluated numerically 
using a L6-poin1 Gaussian quadrature Formula. Eight significant digits were used in the 
calculations and the results were rounded to six significant digits. 

A different method for calculating the values of the integrals corresponding to K=— 40 
is presented here to check t lie values given in table 2. Due to the technique used in preparing 
the table, this will also be a check on the accuracy of the remaining tabular entries. First ob- 
serve that 

/»1q4 i /» oo i 

j[-Ei(-t)]dt= j[-Ei(-t)]dt-O.0OOl+e 



where |«|<0.5(10) _8 , and 






where |^|<1() l0 . [f we use the Maclaurin series for e l , we gel 

From Erdelyi et al. [1953], we find that 

,r(w f l)+Ei(-l)_T(n)-y(n,l)+Ei(-l) 



f 



t n - l [-Ei{-t)]dt-- 



n n 

where T(n,x) and y(n } x) are the incomplete gamma functions, and from Le Caine [1948] we get 



r 



-Ei(-t) 



7 2 + 



dt~- 



6 £ (-1)' 



Thus we find that 



I 



t 2 £=&n\(n+l) 3 

io4 __m:(_f.) 

'' dt~ 0.097843199 



n * 

with an error less than 0.5(10)~ 9 in magnitude. In the same way we may also show that 

*i2 -Ei(-t) 



s: 



dt~ 0.0978432" 



n t 

with an error less than 0.5(10)~ 7 in magnitude. Omitting the details, we also get from the pre- 
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ceding equations that 



s, 



104 pt T 2 



% y(n, 1) 



-Ei(-t)]dt~0M7M3199+j+Ei(-l)[E*(l)-y\-^ ^^-0.0001^0.74509596 



with an error less than 0.5(10) 8 in magnitude. It is seen that six significant digit accuracy is 
indicated for the values of the integrals given in table 2. 



Thanks are due to K. E. Wilkerson for checking the mathematics in this paper, and to Mrs. 
Gail J. Stifel for helping with the computations. 
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